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Abstract 


A method has been described for decomposing 
an optimization problem into a set of subproblems 
and a coordination problem which preserves cou- 
pling between the subproblems. The decomposition 
is achieved by separating the structural element 
optimization subproblems from the assembled struc- 
ture optimization problem. Each element optimiza- 
tion yields the cross-sectional dimensions that 
minimize a cumulative measure of the element 
constraint violations, assuming that the elemental 
forces and stiffness are held constant. The 
assembled structure optimization produces the 
overall mass and stiffness distributions optimized 
for minimum total mass subject to constraints 
which include the cumulative measures of the 
element constraint violations extrapolated 
linearly with respect to the element forces and 
stiffnesses. 


M 

NE 

n 

"(e) 

P 

Q 

STOC 

t(e) 

X 


y 


overbar 


lower and upper bounds, 
respectively 
mass or moment 
number of elements 
length of vector 

number of y-variables in element e 
load in eq. (7); concentrated 
force in numerical example 
elemental force vector 
acronym for subject t£ £onstraints 
number of elemental properties in 
element e 

vector of elemental properties, 
which are design variables at the 
system level 

vector of detailed design variables 
of length n at the subsystem level 
loading case superscript, z = l+NLC 
(number of Joading £ases) 
denotes optimum values 


The method is introduced as a special case of 
a multilevel, multidiscipl inary system optimiza- 
tion and its algorithm is fully described for 
two-level optimization for structures assembled of 
finite elements of arbitrary type. Numerical 
results are given for an example of a framework to 
show that the decomposition method converges and 
yields results comparable to those obtained with- 
out decomposition. It is pointed out that optimi- 
zation by decomposition should reduce the design 
time by allowing groups of engineers, using 
different computers to work concurrently on the 
same large problem. 


Nomenclature 

A cross-sectional area 

C cumulative constraint 

e equality constraint, element 

subscript or superscript 
f functional relation 

F objective function 

g inequality constraint vector of 

length m 

g e elemental inequality constraint 

vector of length m(e) 
g s system inequality constraint 

vector of length m(s) 

K stiffness matrix 

I cross-sectional moment of inertia 

about cent roi dal axis y shown in 
fig. 3 (inset) 


Other notations defined in the text. 


Introduction 

The application of formal optimization 
techniques to the design of large engineering 
structures such as aircraft is presently hindered 
because the number of design variables and 
constraints is so large that the optimization is 
both intractable and costly and can easily satu- 
rate even the most advanced computers available 
today. A remedy is to break the problem into 
several smaller subproblems and a single coordina- 
tion problem, the latter being formulated in a 
manner which preserves the couplings between the 
subproblems. In addition to making the problem 
more tractable, this approach would be compatible 
with the organization of a typical design office 
in which diverse engineering groups work concur- 
rently on different parts of the problem. Such an 
approach would also lend itself to parallel or 
multiple computer processing, thereby shortening 
the design cycle time. 

Several procedures for breaking large struc- 
tural optimization problems into subproblems have 
been proposed in the literature. A typical effort 
is represented by ref. 1 which describes a proce- 
dure consisting of an analysis of the structure 
followed by optimization of each substructure 
while holding invariant the forces acting on it 
from the contiguous substructures. Since the 
optimizations change the stiffnesses of the sub- 
structures, analysis of the assembled structure 
has to be repeated to update the forces acting on 
the substructures for the next round of substruc- 
ture optimizations, and so on, in an iterative 
manner. A similar approach was formulated in 


A'S’S -9C9QS* 



ref. 2. Although computationally efficient thee* 
approaches do not subject the overall stiffness 

algorithm 1 ™ 10 bhe optimization algorithm. That 
tit * th ®refore, cannot be guaranteed to find 
the minimum structural weight because, in qeneral 
a controlled trade-off of the structural material’ 
among the substructures is necessary to find such 

IZTTu i eel “” “ ,s ’*« control,^ 

trade-off, the methods of refs. 1 and 2 are 
basically generalizations of the Fully Stressed 

control "nf h ?h* \ me * hod des igoed to incorporate 
f ?ni to th ! maEena1 distribution among the 

IffeSd n m rof S ?V n aS * emb1ed structure has been 
ottered in ref. 3 for a two-level optimization. 

^ be optimization schemes cited above are all 
rather specialized and would not be suitable for 

wll C o t ’° n to . fflul 1:1 disci pi inary optimization of 
large engineering systems. Recently, ref. 4 

proposed a method for decomposing a large multi - 

«li P ^X t1lrt ^ ti0 " ? roble * 1nt0 a i of 

small subproblems and provided a blueprint for 

methln^Th ° f f h C T5 Ut6r implementation for the 
method. The method decomposes a large problem in 

SeoiaeTbVaT in Each ^bpro P b?em 

dep ’^? d 3 box 1n fl 9- 1 is meant to represent 
a physical subsystem of the total system e a 
airframe or engines in an aircraft/so that'the 
method is entirely general and admits various 
engineering disciplines for analysis of the system 
" d he subsystems. In a particular application 
nHrnnr UCtUreS ’- the dec °mpos i t i on for optimization 
substri!rt C °- nC1 / eS <; Wlth 3 genera1 * multilevel 

analilif so n ?hir re fh* 5 ’ and 7) in stru ctural 
analysis, so that the system acquires a meanina of 

a compete structure, the subsystems bec^ 9 

levei i tUreS> and the subs -ystems of the lowest 
ie y ei, j-j max , correspond to the individual 

finite elements by which the structure is 

ref al 4 Z X‘ln I ?h th1S appllcation > the method of 
rer. g is in the same category as ref. 3 but 

n'rnM™ 1n , th3t n allows 5eve ^ levels of sub- 
problems and in the way the levels of optimization 
are coupled through the concept of optimum 
sensitivity to problem parameters given in ref. 8. 

th „ 1 Exec Jf t ’ on of the method would proceed from 

constraint ^L UPMrd by: W minimizing the 

?ocal d«L • u? ln each sub P rot >lem using its 
local design variables while the higher level 

iaicufitinn r th he1d 3 - “ nstant Parameters, (2) 
calculating the sensitivity derivatives of thi 

subproblem minimum solution to the parameters 
(3) optimizing the system for minimum mass subiect 
to constraints which include the subproblem J 

respect 1 to ^° lat1ons extrapolated linearly with 
respect to the parameters, and (4) reoeatinn 

Sr (1) t,,ro “ 9h (3 > ™ tn is 

.® ecau ?e ref. 4 was intended to only provide 
the blueprint for development of a computer code 
no numerical substantiation of the method was ’ 
included. The method has recently been imple- 
mented for two-level optimization and applied to a 
framework structure as a prelude to proceeding 

onfin.r P J® mentat1on of a 9 en eral multilevel 9 
optimization procedure. The purpose of this paper 

nenor ? escnbe the two-level procedure for the 

nf n ! r h-* case 0f a structu re modeled by an assembly 
of arhitrery type finite elements, and to y 

fi a S c^ at ? ’5 S valid1t -X in a numerical application 
to a simple framework structure which incudes a 


comparison with the results of a conventional 
one-level optimization. invent lonal , 

Two-Level Optimization 

Definitions 

finifl°li the purposes of structural analysis bv a 
sectional S d? nt meth ° d ’ 006 def1nes " cross y 

y = {yi} ; i = i +n (1) 

that can also be organized into NE partitions 

y ~ y e y^ } (2)* 

X t e . collected in a vector X* which is a 

partition e of a vector X for all elements In 
further discussion, the quantities X e are 

referred to as elemental properties Thev ara 
computable as functions of y e : * y 


X e = fe (y e } 


(3a) 


X = { X 1 X® X NE } 


(3b) 

Examples of Xj are: cross-sectional area, 

A, and moment of inertia T a , 

pol.r .oment of , ' J V' 0 ; llT ‘ 

S?hl? 9 “'"I'” “o'Kctent, 0„" ofan 

rWj= PU “- °“ for ,ls«e„t 


M e = f®(X e ) 


(4) 


and stiffness matrix referred to the alnh^i 
coordinate system. Each entry of ‘fhe^SJx K e 


K Pd “ <**> 


(5) 


to. • 

made 'of e1ement * e * could be 


^TT^enotes a column vector written on a single 
nne to save space y,e 
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variable, in which case the material properties 
would be included with the cross-sectional 
dimensions in vector, y e. However, the variable 
material case is outside of the scope of this 
report. 

Although the element mass and stiffness 
appear in eq. (4) and (5) as functions of X e , 
ultimately they are functions of ye through eq. 
(3). Consequently, the elemental properties X e 
and the finite-element cross-sectional dimensions 
y are hierachically related as shown by a Venn 
diagram in fig. 2. The vector y e carries 
information which, for given f e and a set of 

F pq * sufficient to calculate mass and 

stiffness for element e, while the vector X e 
carries the information needed to quantify mass 
and stiffness of the entire structure. 


Proceeding from an element to the assembled 
structure, its stiffness matrix K is generated as: 


single problem at the assembled structure level 
and NE subproblems, one for each finite element, 
at the lower level. The two levels are referred 
X? " s y stem anc * subsystem levels, respectively. 
With respect to a general, multilevel sub- 
structuring scheme shown in fig. l, the two-level 
case corresponds to the upper two levels of the 
figure, with the "substructures" of the second 

HnUe a ele U Ie r n?s. the mean1n9 ° f indlvidual 


Conversion to a two-level optimization scheme 
g 6 ?nto Wlt ^ P art ’ t ’ onin 9 the vector of constraints 


9 = {9 S 


„e 

*9 » • • • 


•9 NE } 


( 10 ) 


where g s contains the constraints on the system 

to h eIrh r ’i and c remainin 9 constraints are local 
to each element. Examples are a nodal point 

fnr P th£ e i"i^ f ° r t !? e f ? rmer and an elen, ent stress 
fnr ns e *u tt6 u‘ Tracing the functional relations 
for g through eqs. (7) and (3b) one obtains: 


K = S(K e ) 
P9 


( 6 ) 


where S symbolizes a procedure for direct sum- 
mation of stiffnesses. Formation and solution of 
the load deflection equations for displacements u 

Ku = P z , z = 1 + NLC ( 7 ) 

where superscript z refers to a loading case, 
yields displacements u and elemental forces 
Q p ’ For element e. Expressed in element 


9 S = f3(X) 


(11) 


? he trace trough eqs. (8), 
( b )> (5), and (3a) leads to: 


g e = f ®(y e .x e ,q e ) 


( 12 ) 


becJlLT" 8 ’ the structural mass * F - In eq. (9a), 


F = f 5(X) 


(13) 


coordinate system. 


Q 


e,z 


K e u 


e,z 


(8) 


where each vector, Q e > z contains r(e) forces 
Q** . The forces in vector Q e » z are 

statically independent and are related to the full 
set of elemental forces by the matrix H e that 
represents the element equilibrium. 

One-Level Optimization 


A conventional, one-level optimization for 
minimum mass can be based on the quantities 
defined by eqs. (l)-(7). Namely, taking y as 
the vector of design variables one has 


min F(y) 

[y } 


(9a) 


subject to constraints (STOC): 
9j(y) < 0; j = Un 


(9b) 


yi < y < y u 


(9c) 


where constraints, g, are imposed on the stati 
behavior variables, such as stresses and displa 
ments, and , y u are side constraints. 


Two-Level Optimization Procedure 

Under a two-level optimization approach, the 
problem defined by eq. (9) is decomposed into a 


vug 

summed. 


maabeb expressed by eq. 


u r c 


a ( e1e ment: ) ^ eve ] *- For conversion to 

* two- 1 eye I optimization, it f s necessary for each 
finite element, e, that the number n(e) of its 
cross^ectional variable 5 be no less than the 
number t(e) of its elemental properties, X e 

t* 


t(e) _< n(e); e * 1 ♦NE 


(14) 


itc Jnnl i^° Ve ! quation b e satisfied in both 
Us equality and inequality parts, or only in its 
equality part, dependent on the type of structural 

pn en 7iai h S i^° Wn by examples in Appendix A. If 

dossi'mI dS 10 1tS ine 9 ua11t y Part, then it is 
possible to carry out an isolated, local operation 
or changing values of the entries in g e by P 

a a wIv U Ihat n ?e the ^ desi9n vaMab1es In ye in such 

constat 1 ?„ ?u’ cons equently, Qe.z reni ain 
on S nai'- * ° ther words > if th e inequality in 

nroni^ L 1S i rue ; there is desi 9n fr eedom to 
proportion the element in a new way. Improved in 

some sense, without affecting the assembled 

«^^t C 3 U^e ^.* 0 • Ut10n • Tra nslated into a formally 
stated optimization problem, that means for y 
element e: 


min C e (g e ) 

(y e } 

* e - F ®(y e ) = 0 


(15a) 

(15b) 
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y? < y e < y e 
l - - u 


(15c) 


Constraints of the problem are the side 
constraints on y as In eq. (9c) and equality 
constraints to enforce Invariance of X e for the 
duration of the solution of eq. (15). Regardless 
of the technique used to satisfy the equality 
constraints In eq. (15b), their presence has the 
effect of reducing the number of free variables in 
y e by t(e), hence, the condition of eq. (14). 

The problem's objective function, C e , is a 
single number that measures the degree of 
constraint violation for all constraints that make 
up vector g e . The quantity C e is known in the 
literature (refs. 4, 9, 10) as a cumulative 
constraint and can be formulated in a number of 
ways. In general, the C e (g e ) should be such 
that: 


>0, if g® > 0; j (l...m(e)); 
(at least one violated) 

C e = < 

<0, if 9® < 0; j = 1 +m(e) ; 

^ (all satisfied) 
and must have continuous derivatives. 


(16) 


Two specific formulations for C e will be 
discussed later (eqs. (23) - ( 25) ) * 

The choice of C e for the objective function 
in the element optimization subproblem, eq. (15), 
is consistent with eq. (4) which for constant 
X e , and for eq. (14) holding, renders F 
unaffected by changes of y e . It means that, as 
similarly proposed in ref. 3, there is no control 
of the objective function at the lower level of 
optimization; the only objective of optimization 
at that level Is to achieve the best possible 
satisfaction of constraints consistent with the 
element forces Q e » z and the elemental properties 
x e . 


System (structure) level The objective 
function and the remaining constraints, g s , are 
controlled at the assembled structure level. The 
single optimization problem to be solved at that 
level is then: 


min F(X) 

{X } 

ST0C gj < 0; j = 1 -hii(s) 


(17a) 

(17b) 


C e < 0; e = 1 +NE 
yf < y e < y e ; e = i + ne 

1 ~ “ u 


(17c) 

(17d) 


variables. Presence of the constraints on C e in 
eq. (17c) and y e in eq. (17d) assures that when 
a solution to the system level optimization 
problem is found, satisfaction of all the local 
constraints will be a part of that solution. 

Coupling between the levels , - The optimiza- 
tion ‘problems in the form given by eq. (15) at the 
subsystem level are coupled to the system level 
problem. It is a two-way coupling: as input, 

each subsystem problem receives the system level 
variables X e and the system level analysis 
results Q e * z (eqs. (12), (15), (16)) and returns 
its optimal C® and y e (eqs. (17c) and (17d)). 
However, in practical implementation the 
constraints on C e and v e cannot remain in the 
form of eq. (17c) and (17d) because their 
evaluation for each new X would require a new 
solution to the subsystem problem in eq. (15)~a 
reoptimization of the elements affected by new X. 
Computationally, that would be prohibitively 
costly in large problems. 

There is a way to bypass these costly 
subsystem reoptimizations in the system level 
optimization. It is available in the concept of 
the sensitivity of the optimum to problem 
parameters and an associated algorithm for 
computation of the derivatives to quantify that 
sensitivity proposed in (ref. 8). Applying that 
concept to the optimization represented by eq. 

(15), one recognizes immediately that the 
optimization constant parameters are X e and 

t 

Q®» z (eq. (12)). Consequently, the optimum 

solution, T5 e and 7 s * of eq. (15) is a function 
of these parameters and has derivatives 
ctfe/dxe dC" e /dQ e > z , d7 e /dX e , and 
t r t 

dy e /dQ^* z , termed optimum sensitivity 
deri vati ves. 

The algorithm of ref. 8 is based on differen- 
tiation of the Lagrange equations that hold at a 
constrained minimum with respect to the problem 
parameters. This leads to a set of linear, 
algebraical equations in which the optimum 
sensitivity derivatives appear as unknown and 
whose matrix of coefficients and the vector of 
free terms include first and second order 
derivatives of the behavior variables with respect 
to the design variables (y e in this application) 
and parameters. Since the computational cost of 
the second derivatives may be significant, it is 
of interest to know that a version of the 
algorithm without the second derivatives of the 
behavior is given in ref. 11. Alternative means 
to reduce the cost of computing these derivatives 
are proposed in refs. 12 and 13. 


Xi < X < X u 

where the objective function: 

F = i M e 
e 


When the optimum sensitivity derivatives are 
(17e) available, they can be used in a Taylor series to 
convert the nonlinear dependence of C e and y e 
on X in eqs. (17c) and (17d) into the linear 
extrapolation approximations: 


depends on X through eq. (4), and the entries of 
the vector X are the system level design 
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(20a) 


9 



ax: 


, e e 
( X t - x to> 


min F ( X) 

{X} 

STOC g* <0 j = 1 *m(s) (20b) 


*1111 

z e t r 


_e e,z 

3Qr /v e e , 

3Q® ,Z 3X® ^ Xt0 ^ 


(19a) 


e 

r = 


+ 1 

t 


if 


3x: 


( x t + V 



3y e e e 
.<• ,x ‘ • x '°> 


• nn 

z e t r 


ar e 

3Qp ,z 


e,z 

3Q r e e 
<X ' ' 


(19b) 


where subscripts a and o denote, respectively, the 
approximate and exact_ values. The above extrapo- 
lation turns TT a and y a into linear functions 

of X only; the dependence on Q e > z is accounted 

r 

for by the chain di f ferentiation in the third term 
in each of the two equations. This chain 
differentiation reflects the dependence of Q e » z 
on X that occurs in redundant structures. In such 
structures, generally speaking, each Q e » z 

r 

depends on each X®, therefore, the summation in 

the chain differentiation spans the entire vector 
X. The derivatives dQ e » z /dX e are in the 
r t 

category of behavior derivatives that are 
routinely available through analytical techniques 
applied to eqs. (7) and (8) (refs. 14,15,16). 

The relations established in eq. (19) will be 
referred to as a linear representation of the 
subsystem. 

In summary, the two optimization levels 
couple through the information flowing between 
them as shown in Table 1. The data passed from 
the system level to the subsystem level carry the 
information defining the X quantities that were 
input into the system level analysis, and the 
output of that analysis in the form of the 
finite-element forces Q e * z . The data 
transmitted in the opposite direction conveys the 
information on the subsystem optima and their 
sensitivity to the data received from the system 
level, for use in the subsystem linear 
representations at that level. 


30 e,Z 

3 ce 9 "r 


z e t r 3 Q e ’ z 3X® (Xt ‘ Xto) 1 ° : (20c) 




= y ♦ I << ♦ 4) 


e = 1 >NE 


e,z 

3 C e 9 Qr 

mi (x* - <.><* (20d) 


2 e t r 3 Q^’ 3X t 


e = 1 -NE 


X 1 < X < x u 

f6(yi>y u ) < x < f 7(yi.y u ) 


(20e) 

(20f) 


x"<X<xJ (20g) 

Two new groups of constraints appear in 
eqs.(20f) and (20g). The constraints of eq. (20f) 
are added to keep the optimization algorithm from 
generating such combinations of X® values that 

cannot be physically implemented at the 
finite-element level (for an example, see eqs. 
(C1)-(C5)). The constraints in eq. (20g) 
introduce additional bounds on X as move limits, 

\ and X M , needed to control the linearization 
1 u 

errors. 

Two-level procedure algorithm and salient 
features .- (he linearization of eos. (17c) and 
TT7d) resembles the local 1 inearization technique 
based on the behavior sensitivity derivatives that 
are known to be very effective in nonlinear 
mathematical programing (refs. 16,17,18). 


System (structure) level problem with 
embedded coupling .- Substituting eg." (191 into 

eq. (17), the system level optimization problem 
becomes: 



in /^? 1 l y, ,- hat bechl ” c l ue could also be used in 
eq. (20b) to linearize g*; whether to exercise 
this option is a problem-dependent decision. As 
it Is the case with any linearization technique 
applied to solve an intrinsically nonlinear 
problem, an iterative procedure has to be 
constructed to allow recovery from the 
linearization errors and the error controlled by 
m0Ve 1 ! fn1tS (eq * ( 2 °9))- In the case 
XlSg steJsT 6 al9 ° rithm C ° nSiStS ° f the 

1* initialize y 

2. compute X (eq. (3)) 

3. analyze assembled structure, obtain Q e > z , 
g , and their derivatives with respect to 
x t ( e qs. (7), (8), and appropriate gradient 

calculation technique) 

4. solve subsystem optimization, eq. (15), for 

each element 1 '* 

5. Calculate optimum sensitivity derivatives for 
the optima found in step 4 

6. Solve the system optimization, eq. (20) 

/. Update X and repeat from step 3 until a 

converged solution is obtained 

In this procedure, optimizations performed in 
' in ?u 4 3re 1 *f rative within themselves and nested 
n the overall iteration spanning steps 3 and 7. 
in further discussion, the latter will be referred 
to as a cycle while the term "iteration" will be 
used in conjunction with step 4. 

a . c Sl 2 C fu the calcu1at ^ons performed in steps 4 
and 5 of the procedure are executed separately for 
each finite element, they can be carried out 
concurrently using distributed computinq 
technology. a 

Information flow between the two levels of 
the procedure is restated in Table 1. Readers 
fami lar with system analysis as formulated in the 
discipline of operations research (ref, 19) will 
recognize the information returned to the system 
level as a particular means to solve the so-called 
system coordination problem. 

In the two-level procedure, the system 
objective function (e.g., structural mass) is 
entirely controlled at the system level bv 
variables X which can be regarded as generalized 
c? ff ,ab that determin e the structure mass 
^ The system objective 

function is not directly included in the subsystem 
optimizations whose only purpose is to achieve the 
best possible satisfaction of the local 

from^hl"! 5 ^ ons 1 1ste ? t with the parameters imposed 
nnon system level . The procedure is entirely 
open to accommodate the designer's judgment as to 

(rlf^PnwJ T i ar ? 6V1Ce of variable linking 
2 be I reely used at both levels to 
keep the number of design variables as small 

r«frai S1 f 6 a ? d ». fo r the same Purpose, one may 
e rain from including all the available elemental 
properties in the set of design variables X (see 
example of a composite panel in Appendix A). 

Equivalence to a one-level optimization 
AcconJTny to eqs. (2u), (lb), and P (15)! the iwo- 

(ref 14? C «nr e ’ Wh6 a conver 9 ed in the Kuhn-Tucker 
(ref. 14) sense, produces a feasible design, just 

as a single-level optimization (eq. (9)) does. 


Moreover if the Kuhn-Tucker conditions are 
satisfied in the subsystem and system level 

that t hL 1 n th ® two : leyel procedure, one can infer 
i % are also satisfied in the y-space in the 
one-level optimization (eq. (9)). Discussion 
the inference is given in Appendix B. In these 
u^ ec ^ s bbe ' ! > tbe two procedures are equivalent. 
However, it does not follow that both will lead to 
the same design point in nonconvex problems havinq 
multiple local minima. In such problems that 
include many practical applications ? the ’soluS ion 
depends on the computational path through the 
design space, and the path taken, in turn depends 
on the algorithm. Since the two procedures are 

rIi^fc hmiCa y d1tferent * a difference in their 
results in nonconvex applications should be 

cXpcCC6u * 

well S’itfrn^ ° f the p [ ocedu,e performance, as 
1tS convergence characteristics and 
overall computational behavior, can only be 
assessed by numerical experiments. Such 

sections"** ^ described in the next two 


Fram ework Structure as a Test Case 
A portal framework shown in fiq. 3 is an 

™tKz J?V hierarchicdl s * stem tb at can be 

optimized for miminum mass under static load 

usinn C L t0 i Stren9th and d1s P |a cement constraints 
using the linear decomposition approach. The 
decomposi ti on is two- level and results directly 
rorn the fact that one can use an engineering beam 

forces fthe na the framework for internal 9 
forces (the end forces on each beam) and 

displacements, assuming that A and I for each beam 

are given but without knowing the detailed cross- 

Si2s°ran s' enS1 ? nS (bl * tl ’**- ) ' TheSe dime "- 

lono as the Pn°S r iZed Separate1 * for each beam as 
long as the end forces in each beam are known and 

sSt f ; x , e v? ,c ?- ,n w™. 

constant A and I of each beam. The correspondence 
of the basic elements of a two-level decomposition 
TablTT t0 the framework exam Ple is given P in * 

and Its One-Level Formula™ 

The framework is composed of three I-beams 
made of the same material and having cross 
sectional dimensions as shown in the inset in 
o’r’fnr Stru ? tural optimization is to be carried 

stL> rfs" 11 " 1 ™" 1 maSS Subject t0 constraints on 
static response induced by two loading cases 1 a 

concentrated force and a coricentrated 9 momenti The 
constraints are Imposed on the framework displace 

l’oa t ifH~ h0 ^ 1 t 0n f a1 translat i° n and rotation at the" 
oaded point of the framework, and on the stresses 

" h eac b — The extreme normal stresses caused 

shoa 6 mo "! ent and axial force, and the extreme 
shear stress due to the transverse force are 
constrained at both ends of each beam to stay 

crit^ca^stresses’of ^ local* 1 buckHng.^he^atter* 16 

s. n ,sts TJSL&r 
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Constraints Include the bounds on the design 
variables. Detailed formulations of all the 
constraints are provided in Appendix C. 

anH i?I'! lyS w 0f the fr fnework for displacements 
and internal forces employs a standard, 
displacement-based, finite-element method 
representing each beam by a single beam element. 
Beam stresses are calculated according to the 
engineering beam bending theory. The critical 
buckling stresses are computed for each part of 
the beam, e.g., a flange, as for an isolated plate 

KtoJSSTrS! 2 b ?r* ry ««<"« routine 

Taking the detailed cross-sectional 
dimensions of each beam as the design variables of 
the problem, the vector y in eq. (2) has 3 
partitions, each containing 6 variables for one 
beam, for a total of 18 design variables: 

y = { { ,t i ,b2» . •} 1, { bittfibg.. .} 2 , 

{bi.tj.b2... } 3 } (21) 


Optimization for minimum mass can be replaced 
with optimization for minimum material volume 
because of the material homogeneity. Denoting the 
beam length by 1-j , the objective function 
becomes a function of cross-sectional area A: 

, F = M “ Z Ajli (22) 

i 

S Mnn P ^h em “ n J be . S0,ved as a one-level optimize- 
tion with n=18 design variables in a conventional 
tormulation such as given in eq. (9). 

Two-Level Formulation 

Under the two-level approach, the framework 
system) is considered decomposed into three beams 
(subsystems, finite elements) under the action of 
the beam-end forces shown in fig. 3 (inset), if 
the decomposed framework were superimposed on the 
general, multilevel decomposition scheme shown in 
ng. l, the assembled framework would fall in the 
J! n ii re ? t ™'-ture" bo * at level 1 and each beam 
would coincide with a substructure at level 2 
although in this case "substructure" simplifies to 
a single finite element. For the purposes of the 
framework analysis, each beam's stiffness and mass 
properties are determined by two elemental 
properties: cross-sectional area A and moment of 

inertia I. These quantities become the system 
level design variables in vector X, eq. (3b) as 
indicated in Table 2. The framework displacement 
constraints are in the g s category, while the 
beam stress constraints are included as g e in 

Since the condit1 °n in eq. (14) holds 
(t(e)-2<n(e)=6) (see Table 2) for each beam, the 
original problem can be solved decomposed into 
three subsystem problems of six design variables 
y each and a system problem of six design 
variables X, according to eq. (15) and eq. (20), 
respectively. M v ' 

Solution of eq. (20) coupled with eq. (15) 
requires a specific form of the cumulative 
Constraint C« in eq. (20c). Two different 
formulations for the cumulative constraint were 
studied in the test case. One of them is 


suggested by the concept of a quadratic exterior 
penalty function: 

c ® = I. «9j>) 2 (23a) 


V - 


9j . if 9j > 0 

0.0, if g < 0 


For continuous g,- functions, this formulation 
has continuous first derivatives which permits use 
of gradient-dependent algorithms in solution of 
eq. (3). Another function that can be used for 
the same purpose is a function proposed in ref 22 

1 ?. to a PP™ximate the maximum 
constraint. The function, henceforth referred to 
by the acronym KS, has the form: 


i m 

K$(g . ) = — In ( l exp 
J J = 1 


( p9 )) 
J 


coJsfraiS 6 PPOperty of f °1 lowing the maximum 
MAX( 9j ) < KS < MAX(g ) +1 l n ( m ) (2E 


with a tolerance that depends on the constant o 
supplied by the user. The KS function, like the 
quadratic exterior penalty function, has 
continuous derivatives; in addition it performs as 
an extended penalty function because it is defined 
throughout the infeasible as well as the feasible 
domains. 

The seven-step iterative procedure for 
solution of the two-level optimization listed 
previously has been implemented for the framework 
test case in the manner described in Table 3. 


Numerical Results 

Two-level structural optimization by linear 
decomposition is demonstrated for the framework 
example. The procedure given in Table 3 is 
implemented in a Fortran main program that calls a 
finite-element analysis subroutine based on a 
stiffness method and an optimization subroutine 

&iM m ref - 23) that ^p 10 ^ a usable - 

I IIImc 7cn reCt1on techni 9 ue - Results obtained on 
a KKiMt 750 computer include benchmark data for a 

conventional, one-level optimization and the 
two-level optimization data. The detailed 
formulation of the constraint functions, including 
equality constraints on Aj and I,* used in all 9 
the numerical tests are given in Appendix C. 

Res ults for a Conventional, One-Level Approach 

Several variants of the optimization and 
different starting points in both the 
feasible and infeasible domains were used to 

was a chosJn re * u1ts - In variant 1, which 

was chosen to be the reference technique, all 

constraTnts were kept separate, while in variant 2 
the constraints, except side constraints, were 
collected in a cumulative constraint. A 
piecewise-1 inear procedure combined with the 
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cumulative constraint (in the form given in eq. 
(24)) was carried out in variant 3 using move 
limits of 15 percent. The purpose of including 
variants 2 and 3 in the benchmark testing was to 
determine to what extent the optimization results 
were influenced by the use of a cumulative 
constraint and a piecewise-linear procedure, which 
are both embedded in the two-level optimization. 

It turned out that the three variants and 
different starting points generate designs having 
masses which fell within 5 percent of the variant 
1 result. However, there was as much as a 300 
percent difference in some design variables. 

The dependence of the optimum on the starting 
point and search path indicates that the problem 
is nonconvex and has local optima, and that a weak 
functional relationship exists between the 
objective function and constraints and at least 
some of the design variables (a "shallow" 
optimum). In this particular example, the scatter 
of the local minima is bounded by two extreme 
cases both obtained by variant 1 starting from a 
feasible design and from an infeasible design near 
minimum gage. The two local minima differ 
significantly in their distributions of A and I 
shown in fig. 4 and also their local dimensions 
given in Table 4. The optimum design shown in 
fig. 4a transmits the load to the ground support 
primarily through flexural stiffness of the right- 
hand side vertical beam. In the optimum design 
depicted in fig. 4b, the load is transmitted 
primarily through flexural and extensional 
stiffness of the horizontal beam to the left-hand 
side vertical beam which is relatively stiff in 
bending due to its shorter length. Despite the 
differences reaching 58 percent in I, the 
structural masses differ by only 6 percent, and 
all constraints are satisfied. The deliberate 
exclusion of the displacements out of the plane of 
the framework eliminated beam torsional -bending 
buckling and resulted, predictably, in the 
unusually large depth-to-flange-width ratio 
indicated in Table 4. However, since these 
proportions have no bearing on the purpose of the 
numerical verification of the method, the two 
bounding cases of variant 1 were selected as 
acceptable benchmark results. 

Results for the Two-Level Procedure 

Numerical studies with the two-level 
optimization were carried out for two formulations 
of the subsystem cumulative constraint: the 

quadratic exterior penalty function defined in 
eq. (23) and the KS function defined in eq. (24). 

Cumulative constraint in form of an exter i or 
penalty function .- Optimization with thp 
exterior penalty function formulation (eq. (23)) 
consistently yielded designs about 15 percent 
heavier than the reference result. It was clear 
from examination of the history of iterations that 
the discrepancy was caused by the loss of the 
gradient information near and in the feasible 
domain where the exterior penalty function and its 
derivatives vanish. The gradient information was 
needed by the optimization program (which is based 
on the usable-feasible directions method) and the 
iteration history showed that as long as that 
information was available the optimization program 
successfully guided the growth of the structure 


from the near-minimum-gage starting point toward 
the feasible domain. However, it was unable to 
reduce the structural mass while maintaining a 
feasible design once the feasible design space was 
entered. It was concluded that the quadratic 
exterior penalty function is not a proper choice 
for a cumulative constraint to be used in the 
context of the proposed multilevel optimization 
method and attention was then directed to the use 
of the KS function (eq. (24)), for the cumulative 
constraint. 

con straint in form of KS 
functi on .- Sat i sTa ct ory results were obtained 
using the KS function. They are collected in 
Table 4 and show the method’s ability to generate 
designs comparable to the benchmark design when 
starting from the same point, either feasible or 
infeasible. In fact, the objective functions 
obtained by means of the two-level optimization 
started from infeasible and feasible points 
exceeded the benchmark value of the objective 
function by only 1.6 percent and 1.9 percent, 
respecti vely. Another way to assess the 
effectiveness of new optimization method relative 
to the reference method is through comparing the 
changes they generate in the objective function. 
Denoting by V the ratio of the final to initial 
values of the objective function, one may take 
Il-rl as a measure of the change relative to the 
initial value. Table 4 shows that the two-level 
optimization overestimated the above measure of 
change by 2.2 percent when starting from an 
infeasible point and underestimating it by 95 
percent when starting from a feasible point. To 
make this comparison meaningful, the starting 
points were deliberately chosen to make the ratios 
"r" significantly different from unity, 3.43 and 
0.329 for the infeasible and feasible starting 
points, respecti vely. 

More significant differences were recorded 
among the individual design variables, at both 
local and system levels, in the optimal designs 
corresponding to different initial points. As one 
may see in Table 4, these differences reach 350 
percent. They confirm the distinctly nonconvex 
and “shallow" optimum nature of this particular 
example problem. Since these are consistent with 
similar differences observed in the one-level 
optimization, they were not introduced by the 
two-level method itself. 

Converg ence hist ory.- Convergence of 
optimization is illustrated by history plots in 
fig. 5 for the one-level optimization and in 
figs. 6 through 12 for the two-level procedure. 
Abscissas of these plots refer to iterations of 
the all-in-one optimization and cycles of the 
two-level optimization. One iteration is one of 
many consecutive executions of the usable-feasible 
directions algorithm (for more precise definition 
see description of program CONMIN in ref. 23) 
while one cycle is one execution of steps 4 
through 11 of the procedure described in Table 3. 

The convergence was found to be similar in 
character to that of the conventional, one-level 
method as illustrated by comparison of figure 6 
with figure 5 for an infeasible starting point. 

In both figures, the objective function rises, 
overshoots the optimal level and then returns to 
it asymptotically. Constraints are reduced to 
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feasible (negative) values as illustrated by an 
individual stress constraint in fig. 5 and the 
cumulative constraint of beam 1 in fig. 6. (The 
factor of 5 in fig. 6 is for scale uniformity.) 
The jagged character of the plots is inherent in 
any algorithm based on the usable-feasible 
directions method and is more pronounced in 
two-level than in one-level optimization for 
reasons to be discussed later. 

Objective function history starting with a 
feasible design which is also shown in fig. 6 
indicates convergence after 18 cycles as compared 
to about 28 cycles for an infeasible starting 
point. An example of the behavior of the system 
level variables is illustrated by the histories of 
the cross-sectional areas of the beams in fig. 7a 
and 7b. It is apparent from comparison of fig. 7b 
with fig. 6 that the design variables require a 
few more cycles to converge than the objective 
function. To complete the illustration of the 
system level optimization history, the framework 
displacement constraints are plotted in fig. 8. 

At the subsystem level, the optimization 
history graphs are given in figs. 9 through 12 for 
beam 1. The behavior in the other two beams is 
similar. Graphs for variables bj, ti, b 2 , 
and t 2 (see fig. 3) are shown in figs. 9a and 
9b. They repeat the familiar pattern of rise- 
overshoot-descend for infeasible starting designs 
and a smoother asymptotic descent for the feasible 
starting designs. The variation of the beam 
objective function (cumulative constraint) is 
shown in fig. 10 which depicts clearly the 
elimination of local constraint violations when 
starting from an infeasible design, and a 
reduction of their oversatisfaction as the beam is 
being slimmed down in the optimization that begins 
from a feasible design. One of the subsystem 
level constraints that makes up that objective 
function is local buckling of the beam flange; its 
history plots in fig. 11 correspond to those in 
fig. 10. Each graph plotted in fig. 11 shows the 
value of constraint at the end of subsystem level 
optimization (step 6, Table 3) that was carried 
out for beam 1 in each cycle. To illustrate the 
character of the constraint changes during the 
subsystem optimization, the constraint history for 
cycle 8 for the infeasible design start case is 
plotted in fig. 12. 

Accuracy of linear extrapolation .- The 
multilevel optimization approach is predicated on 
the accuracy of the linear extrapolations based on 
the optimum sensitivity derivatives. Therefore, 
it is interesting to see how the cumulative 
constraint predicted by the linear extrapolation 
at the end of one cycle (Table 3, steps 10 and 11) 
compares with the result of full analysis carried 
out at the beginning of the next cycle (Table 3, 
step 4). Such a. comparison is displayed in fig. 

13 and shows that the prediction error eventually 
vanishes for a sufficient number of cycles thus 
permitting the procedure to converge. The graphs 
show that before the convergence is reached the 
linear extrapolation consistently underpredicts 
the cumulative constraint value when proceeding 
from an infeasible design starting point and' 
overpredicts it when the start is made from a 
feasible design point. 


Characteri st i cally , the relative error is 
larger when the optimization is started from an 
infeasible design, apparently because the 
procedure then goes through a number of changes in 
the membership of the active constraint set that 
comprises constraints defined in eq. (15c) for 
each beam. Consistent with observations reported 
in refs. 8 and 24, these changes tend to degrade 
the accuracy of the optimum sensitivity 
derivatives as the behavior predictors. The 
larger prediction errors apparently cause the 
history plots to be somewhat more jagged in all 
cases for which the optimization is started from 
an infeasible design rather than from a feasible 
one. They also slow down the convergence, so that 
a larger number of cycles is required when 
starting from an infeasible design. However, the 
procedure exhibits a good ability to recover from 
occasional large prediction errors, e.g., cycle 8 
in fig. 13 (infeasible design start) and cycles 11 
and 13 in fig, 13 (feasible design start), and 
gets back on track with remarkable robustness. 

Computational cost .- Since the purpose of the 
reported work was a demonstration of a concept, no 
attempt was made to refine either the reference, 
one-level procedure or the two-level procedure for 
maximum computational efficiency, especially since 
the framework example is much too small to 
demonstrate efficiency of any optimization or 
analysis method. Nevertheless, one may obseve 
that the two-level optimization converges in a 
number of cycles about equal to the number of 
iterations in the one-level optimization, with the 
numerical workload in a cycle being less than in 
an iteration. While the precise workload 
difference depends on the algorithmic and 
implementation details, the major difference stems 
from having to calculate a number of gradient 
vectors equal to the number of design variables 
which at the system level of the two-level 
optimization is smaller than in the one-level 
optimization. Specifically, in the framework 
example the number is reduced from 18 to 6. 
Computational cost savings resulting from that 
reduction are partly offset by the cost of the 
subsystem level optimizations and the associated 
sensitivity analysis. Approximating computational 
cost by the CPU time and setting the time for one 
cycle at 100 percent, the times spent at the 
system and subsystem levels within one cycle were 
28 percent and 72 percent, respecti vely. Thus, 
the total time share for each of the three beams 
(subsystems) was 24 percent, a nearly uniform 
distribution of time among the system and each of 
the subsystems. These relative cost values are, 
of course, strongly problem dependent as is the* 
total cost of execution of the entire procedure. 

That cost is expected to become smaller relatively 
to the cost of a one-level procedure as the system 
grows in terms of the number of elastic degrees of 
freedom and number of structural components 
(subsystems). 

Concluding Remarks 

A method has been described for decomposing 
an optimization problem into a set of subproblems 
and a single coordination problem which preserves 
the coupling between the subproblems. The 
resulting procedure is iterative and calls for 
repetitive analysis of the assembled structure, 
optimization of the individual components as 
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subproblems, followed by optimization of the 
assembled structure in which the component optimum 
solutions are extrapolated linearly using their 
optimum sensitivity derivatives with respect to 
the system level design variables and internal 
forces. The subproblems are organized 
hierarchically into two levels. The variables at 
the lowest level (the subsystem level) are 
physical cross-sectional dimensions; the variables 
at the highest level (the system level) are 
quantities that govern the stiffness and mass 
distribution among the finite elements of the 
structure. The overall objective function (such 
as structural weight) is controlled at the system 
level. Optimization at that level influences the 
level below by means of changing the mass and 
stiffness distributions, and the associated 
distribution of internal forces. 

The method is demonstrated using a portal 
framework as an example of a two-level structure 
in which the system-level variables are the 
cross-sectional areas and moments of inertia of 
the beams, and the subsystem-level variables are 
the beam detailed cross-sectional dimensions. 
Verification of the method by comparison with the 
results obtained by a conventional, one-level 
optimization show the validity and effectiveness 
of the proposed approach. 

Satisfactory testing of the two-level 
approach is a stepping stone for implementation of 
a multilevel structural optimization procedure. 

That implementation is seen as a stage in 
development of a multilevel optimization for 
multidisciplinary engineering systems whose goal 
is to allow groups of engineers using distributed 
computing technology to work concurrently on 
various parts of the problem, thereby reducing the 
real time of the system design. 
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Appendix A 

Discussion of the Problem Dimensionality with 
Examples 

A stiffened panel of the type frequently used 
in aircraft wing covers is an example of an 
element satisfying eq. (14). Assuming the panel 
cross-section to be as in fig. Al, the number of 
cross-sectional design variables is n(e)=8, if 
stringers are laid in one direction only. To 
represent the panel as an orthotropic membrane 
finite element in structural analysis and in a 
minimum mass optimization, the stiffnesses Axx, 

Ayy, Axy, and mass M are needed which, for an 
isotropic material, are uniquely defined by only 
two elemental properties: TS-skin thickness, and 

TR-equi valent ("smeared") thickness of the 
stringers to be used as the system level design 
variables. Thus, t(e)=2 < n(e)=8. 

If the panel were made of several layers of a 
composite material, laid in, say, four different 
orientation angles then, including the angles 
withassociated thicknesses in the set of the 
cross-sectional design variables, n(e)=2*4=8. The 
panel mass and membrane stiffness definition 
requires, as above, total mass M, and 
orthotropic stiffness coefficients Axx, Ayy, and 
Axy. The total mass is defined by total thickness 
T, but, in contrast to the stiffened panel, the 
stiffness coefficients depend on all cross- 
sectional design variables in an arithmetically 
complex way which does not make it practical 
expressing them in a closed analytical form by a 


subset of the elemental design variables. 
Therefore, the elemental properties to be used as 
the system level design variables are: T, Axx, 

Ayy, and Axy. Hence, t(e)=4 < n(e)=8. 

For an opposite example, consider a beam of a 
solid circular cross-section. Its diameter D is 
the only cross-sectional design variable, n(e)=l. 
On the other hand, for a general case of 
structural analysis and minimum mass optimization, 
the beam must be represented by its cross-section 
constants A, lx, Iy, J, and mass M. Since all 
these quantities are uniquely defined by the beam 
diameter D, that diameter is the only elemental 
property that can be used as a system design 
variable. Thus, in this case t(e)=n(e)*l. 

In closing, let us return to the composite 
panel example to make an important point. The fact 
that there are four elemental properties available 
does not require using all of them as system level 
design variables. Using them all means full 
control over the element mass and stiffness. It 
is easy to imagine applications where such full 
control is unnecessary and giving up a part of it 
for the benefit of reducing the number of system 
level design variables may be a reasonable 
choice. For example, one may choose to use only 
M and Axx as design variables while allowing 
Ayy and Axy to be "passively" computed from the 
panel cross-sectional design variables. 

Appendix B 

Equivalence of the One-Level and Two-Level 
Optimizations with Respect to" the Kuhn-Tucker 
Optimality Conditions 

Suppose that the two-level optimization 
procedure converged to a design in which the 
Kuhn-Tucker (K-T) optimality criteria are 
satisfied at the system level and for each element 
at the subsystem level. A natural question to 
ask at that point is whether it would be possible 
to obtain a better design, i.e., lower objective 
function value, without violation of the 
constraints by continuing the optimization in the 
y-space in a conventional one-level way. In other 
words, the question is whether there is anything 
in the two-level approach that would make it stop 
short of the design that could be generated by a 
one-level procedure. Examination of the K-T 
conditions at the design point to which the 
two-level procedure has converged provides the 
answer. For the two-level procedure, the K-T 
conditions corresponding to the system level 
eq. (20) are: 



e 

C = 0; e = 1 + NE (Blc) 
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X > 0 
J- 

X e ' > 0 


Simultaneously, at the subsystem level the 

with^he’f 10 " 5 /^ each e1ement a re consistent 
with the formulation in eq, ( 15 ): 


: e 3g t 

7 * ‘« sjf * ? 


t b 3y i + b b ^7 * 0; 1 = 1 * n(e) 


0.015 rd, respectively. In addition, there are 
constraints on the system level design variables 
to guard against occurrence of physically 
impossible combinations of these variables (e a 
moment of inertia disproportionately large with*’ 
cros ?- s ® ct1ona l area). These 9 
(20f) 'are: COnstra1nts ’ correspond to eq. 

3 

Ij m1n * bl min t lmin A ^min^min 


9 t * x t " f 4 ( y6 ) = 0; t = 1 + t(e) 


9k - 0 


t3 ,nin h min 


x b > 0 


where g 5 denotes active side constraints on 
variables y . 

awavf^fi^ i Bla) means that no further "wvement 
away from the design point is possible in the 

fhrnrtnh /D^ h \ Ut )fi 0,at,n 9 constraints in eq. (Bib) 
through (Ble). Since eq. (B2b) can be interpreted 
as linking t(e) of the n(e variables P a d 

t e n i^L° the . va ^ ab1 « xe, it follows that 
£n! 1 9ked variables y can not be changed also 
However, the remaining n(e)-t(e) variables y in 

S ildirTt dl fK f an n0t be chan 9 ed because^eq. 

B2) indicates that no further reduction of the 

Violin nn V %K COnStra i nt * <S P°“ible without 
JbSJi 9 riS« constra ! nts in eq. (B2b) through 

changed J!i nA 11006 ° f the y variab1es can be 
ganged and no movement away from the point is 

t*l iVrl 10 ^ deSi9n $ P ace which is exactly 
the same conclusion that would result from * 

Point f foJ’a n on f 1 th6 f‘ T cond i tions at the same 
entirelv in ; leve ' optimization conducted 
that lfth 2 J hat , Spa f e * Hence > one may assert 
Sf Ltisfarf °‘ TfJ procedure terminates because 
ve °? ° f ! he . K - T conditions at both 

point -in the K v T c COOd ^ ,0nS at the corresponding 

satisfied also" Pa Th 0 f °f 3 one ’ leve1 Procedure are 
satisried a so. Therefore, the termination point 

?L?t tW °; ,eVel Procedure must be at least a ? 
local ^constrained minimum. 

noted in the body of the paper, this 
assertion does not mean, however, that the two 
procedures will arrive at the same result 
evenwhen started from the same point if the 
optimization problem is nonconvex. 

Appendix C 

Detailed Formulation of the Framework Test Case 
This appendix provides details of the 

systein d leve(°and ’ at ' t h f< V he framework ab the 
system level and at the subsystem (individual 

« L*!L 

35 :*S 1ST. 


bl min tl inin^ min " — — ^ ) + 

2 


+ b2 min t2 min^ t lmin +hm in + “ 


^2 + 


hm i nt 2 mi n (t n + h|ni n "^i n ) 


v . = bl min tl rnin 

y mln - — + bp . t? . (ti 

2 ‘min *min' Imin 


+ b min + 


x h . f- fa. ' f, min 
m1n ^min^min * — “ )/ 

^min^min 

+ b2 min t2 min + ^n^in) (Cla 

^max " * s ab ° ve » replace subscript "min" with 
max • 

(C2) 

Ai min = bl min t lmin +b 2min t 2min + 

+ ^n^min (C3) 

A ’max ° ab ° ve ’ replace subscript "min" with 
ma x • „ 

(C4) 

! ’max "• i ^B i ^max (C5) 


o 3 

Ii A = bl| " <ntl min + ^min^min + 


t3 m <_h 3 
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s = 2 (AjbjCi + A2B2C2 + A3B3C3) 


+ l>2 m1n t2 m1n^ tl m1n + j - y) 2 + 

* "W'%. * | - >> 2 <«*> 

in which 


[ 


bl min tl mir 


12 


b2 mi n t2 min^ tl min +h 


t2 min 

+ - ) + ^3 m ^ n (ti m1n + *Ly 

^ b l mi n b ^i n +b2 mi n t2 mi n +ht 3 mi n ^ 


(C5b) 


and 


h - (A-bl m j n ti m ^ n -b2 m j n t2 m i n )/ 
t3 min 


(C5c) 


and 


Ii B = [ 1+2 (MOVE LIMIT)] I y . 


(C5d) 


2 2 2 
T =* A^Bj + A 2 B 2 + A 3 B 3 - I m 

I m ■ I - Ajt|/12 - A 2 ti/12 

Ai = b2t2 
A 2 » 

A 3 » A - Aj - Ag 

B 0 = (Ajt 2 /2 + A 2 t 2 + Ajtg + 
+ Agt^J/A 

Bi - B 0 - t 2 /2 
B 2 = B 0 - t 2 - t i/2 
B3 = B 0 - t 2 

C X = (A 3 /2 + A 2 )/A 
C 2 = Cl - 1.0 

c 3 = c x - 0.5 


} 


(C9a) 


) (C9b) 






J 



(C9c) 


(C9d) 


Constraints in eqs. (Cl), (C2), (C3), and 
(C4) relate the upper and lower limits of the 
system level design variables to those of the 
local design variables* The constraints given in 
eqs. (C5a) to (C5d) keep the cross-sectional 
moment of inertia commensurate with the 
cross-sectional area. 

At the subsystem (individual beam) level, the 
constraints are prescribed for an I-shaped 
cross-section as shown in fig. 3 which also 
defines the dimensions used as local (subsystem) 
design variables. Since there are two equality 
constraints: 

ej = Aj/AS-l = 0 (C6) 

e 2 3 li/ 1 *-! = 0 (C7) 

one can eliminate two local design variables by 
expressing them in terms of the remaining, 
independent variables. The appropriate relations 
obtained from eq. (C6) and eq. (C7) are: 

h » (-s + /s 2 - 4RT)/2R (C8) 

t 3 = A 3 /h (C9) 

where 


R = A 3 /I 2 + A^i + A 2 c| + A3C3 (C 9a) 


and the remaining problem is constrained only by 
the side constraints and inequality constraints. 
The expedient of using a closed-form solution for 
the equality constraints (which was used also in 
ref. 25) does not detract from generality because 
any established technique for handling equality 
constraints could have been substituted for the 
closed-form solution if such a solution were 
unobtainable. 

The side constraints are minimum gages and 
upper bounds as follows: 

b m in 3 10* bmax = 100.0, stc. 

The inequality constraints are imposed to 
prevent overstress and buckling. The overstress 
constraints are written in a standard form: 

9j " ^ - 1 ( C10 > 


where the material allowable a a s 20,000 
N/cm 2 and o is a normal stress due to axial 
force combined with a bending moment and computed 
by means of a textbook formula N/A + Mc/I. 
Equation (CIO) is evaluated at four points, top 
and bottom of the cross-section at both ends of 
the beam and, thus, represents four constraints. 
Similarly, a constraint on the maximum shear 
stress is: 

9j - h - 1 (Cll) 

a 
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al I 0Wab ! e *a * 11.600 N/cm2 and 
t - VQ/It. Equation (Cll) Is evaluated at the 

a„T^ S ?, Cti0n Centroid at both ends of the beam 
and so it represents two constraints. 

ct;,M?^ lin ;! consbraints guard the beam flange 
stability under action of normal compressive 9 

2hea? S str.« rrT Web Stab111t * ^r action of 
shear stress. It Is assumed for buckllnq 

evaluation purposes that each flange is a plate of 
I(d.! 1,1 ? ,dth ? supported .long 

sn- rjX^ris'rpfr;,^- L 
s,des - "» sti " a " a 


9j = |° |/°ab - 1 
9j - | T |/ T ab - 1 


where <j and r are computed as for the 
overstress constraints and the critical values of 
stress are obtained from the familiar formulae: 


°ab = 0.4 (0.904) E _ 
b 


/ 1 

Tab = 5*5 (0.904)E _ 

i b 


hn«-h I!! e buck li n 9 constraints are evaluated for 

KEh e t n h d e sTt r he a tlr r * 


Table 1 Interlevel flow of Information 
System level (assembled structure) 


elemental properties 
and forces 

X e ; Q 6 ’ 2 


elemental optimum 
solutions and 
their derivatives 

V; y* 

ac e /3x 

3^/8 Q e,z 
r a 


a P/ax 

ay e /3Q e * 2 

r 


Subsystem level (Individual elements) 


Table 2 Correspondence of the quantities in the 
framework example to the generic quantities 


GENERIC 


FRAMEWORK 


y = {•••{y®}....} y = { ....{b 1 ,t 1 ,b 2 ,t 2 ,h,t 3 } ...} 


x * = f e (y e ) jA | e 

j " f ® l b I* t l. b 2.t2.h.t 3 } 

e = 1, 2, 3. 


tti 

NE x t(e) 

1 (e) 

r(e) 

z 


{ N, M, T }e,z 

constraints on the loaded node 
horizontal translation and 
rotation due to P and M 

beam stress and local buckllnq 
constraints s 

beam mass M = f(A 1 ,A 2 ,A 3 ) 

ri ’ 

3 

18 

6 

2 

6 

3 

3 x NLC =6 
2 
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Table 3 Two-level procedure implemented for the framework test case 

This table contains a step-by-step summary of the two-level optimization 
procedure for the framework. 


System level - whole framework 

Component (subsystem) level - 
each separate beam 

1. Define loads. 

2. Define displacement constraints. 

1. Initialize detailed dimensions. 

2. Compute A, I for each beam. 

3. Move to the system level. 

4. Analyze the framework to compute 
Its displacements and the end 
forces (N,M,T, fig. 2) on each 
beam. Compute derivatives of 

- these quantities with respect to 
the A, I of each beam. 

5. Move to the component level. 



6. For beam 1, hold constant the end 
forces N, M, T and the values 
A and I. 


Analyze the beam to evaluate its 
constraints, such as stress and 
local buckling. 

Form a single measure of the 
constraint violation using, for 
example, an exterior penalty 
function. 

Optimize 6 cross-sectional 
dimensions as subsystem design 
variables to minimize the measure 
of constraint violation as an 
objective function subject to 
minimum gage and other side 
constraints, Including equality 
constraints on A and I. 

The equality constraints assure 
that the beams A and 1 
computed from the cross-sectional 
dimensions are equal to those 
prescribed at the system level. 

/. For optimized beam, compute de- 
rivatives of the minimized meas- 
of the constraint violation and 
the subsystem design variables 
with respect to the constants: 

N, M, T and A, I. 
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Table 3 Concluded. 



8. Repeat steps 6 and 7 for beams 
2 and 3. 

9. Move to the system level. 

10. Approximate the minimized 


measures of the constraint viola- 


tion in each beam as linear func- 


tions of 6 quantities A and 


I by a linear Taylor expansion. 


using the derivatives computed 


in step 7. In this expansion 


each of the end forces N, M, T 


is also approximated as a linear 


function of all 6 quantities A, 


and I using derivatives com- 


puted in step 4. 


-11. Optimize 6 system level variables 


A, 1, to minimize structural 


mass subject to: 


(a) framework displacement con- 


straints approximated as func- 


tions of A and I by 


linear Taylor expansion using 


derivatives computed in 


step 4. 


(b)constraints requiring that the 


minimized measure of the con- 


straints violation in each 


beam be reduced by a pre- 


determined decrement. 


(c)move limits on the variables 


A and 1 to protect 


accuracy of the linear Taylor 


expansions and to account for 


side constraints of the sub- 


system design variables. The 


latter are approximated as 


functions of A and 1 by a 


linear Taylor expansion using 


derivatives computed in 


step 7. 


(d)side constaints on A and I. 


12. Go back to step 4 with the system 


level design variables A and 


I obtained in step 11, and the 


corresponding approximate sub- 


system deign variables estimated 


by a linear extrapolation as 


described in step lie. 


Terminate when: 


(a)the framework displacements 


are within constraints. 


(b)the minimized measure of con- 


straint violation for each 


beam is reduced to at least 


zero. 


(c)no further reduction of the 


framework mass appears 


possible. 
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Table 4 Comparison of Optimization Results 


LEVEL 1 
LEVEL 2 

LEVEL j 

LEVEL 

*max 



INITIAL 

VALUES 

FINAL 
RESULTS 
ONE LEVEL 
OPTIMI- 
ZATION 

FINAL 

RESULTS 

TWO-LEVEL 

OPTIMI- 

ZATION 

INITIAL 

VALUES 

FINAL 
RESULTS 
ONE LEVEL 
OPTIMI- 
ZATION 

FINAL 

RESULTS 

TWO-LEVEL 

OPTIMI- 

ZATION 

Obj. 

26,469 cm 3 

90,682 cm 3 

92,090 cm 3 

275,000 cm 3 

90,592 cm 3 

92,330 cm 

Beam 1 







»i 

11.0 

10.0 

10.3 

30.0 

13.0 

10.3 

ti 

0.275 

0.491 

0.571 

1.0 

0.450 

0.569 

h 

22.0 

78.1 

73.9 

50.0 

74.9 

74.0 

t 3 

0.275 

0.517 

0.518 

1.0 

0.497 

0.519 

b 2 

5.5 

8.08 

5.08 

30.0 

12.1 

5.13 

t 2 

0.275 

0.511 

1.18 

1.0 

0.487 

1.16 

Beam 2 







bi 

11.0 

10.3 

10.3 

30.0 

11.4 

10.7 

t \ 

0.275 

0.439 

0.476 

1.0 

0.404 

0.451 

h 

22.0 

95.7 

89.4 

50.0 

89.9 

90.1 

t 3 

0.275 

0.421 

0.414 

1.0 

0.397 

0.417 

b 2 

5.5 

8.46 

5.14 

30.0 

10.7 

5.12 

t 2 

0.275 

0.539 

0.984 

1.0 

0.435 

0.960 

Beam 3 







bi 

5.5 

5.0 

5.03 

30.0 

7.50 

4.98 


0.275 

0.267 

0.253 

1.0 

0.268 

0.253 

h 

22.0 

47.8 

59.1 

50.0 

61.9 

59.0 

*3 

0.275 

0.25 

0.251 

1.0 

0.25 

0.25 

b 2 

11.0 

10.0 

10.0 

30.0 

10.0 

10.0 

t 2 

0.275 

0.332 

0.400 

1.0 

0.369 

0.393 




NOTES: See Figure 3 for dimension definitions. 

All beam dimensions are in centimeters. 



Fig. 1 Multilevel substructuring. 
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SCALE FOR 
FEASIBLE 


Fig. 3 A portal framework. Fig. 6 History of objective function (mass) for 

the entire structure and cumulative 
constraint for beam 1 in two-level 
optimization starting from an infeasible 
design, and the objective function (mass) 
for the entire structure starting from a 
feasible design. 


INFEASIBLE STARTING 

POINT FEASIBLE STARTING POINT 



Fig. 4 Distribution of A and I for two 
optimal designs of the framework. 


DESIGN 

VARIABLES 



Fig. 7 History of selected system level design 
variables, (a) infeasible design start. 



Fig. 5 History of objective function (mass) and 
one of the flange buckling constraints 
in the benchmark, one-level optimization, 
starting from infeasible design. 
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Fig. 8 History of one of the system level 
constraints— horizontal displacement 
due to force P. (a) Infeasible design 
start. - (b) feasible design start. 


Fig. 10 History of the cumulative constraint 
used as objective function in optimi- 
zation of beam 1, (a) Infeasible design 
start. - (b) feasible design start. 



Fig. 9 History of local design variables for 
beam 1, (a) infeasible design start. 



CYCLfS 


Fig. 11 History of one of the local constraints 
(flange buckling In beam 1), (a) 
infeasible design start. - (b) feasible 
design start. 



Fig. 9 History of local design variables for 
beam 1, (b) feasible design start. 



ITERATIONS 


Fig. 12 Detailed history of the constraint from 
fig. 11(a) over the iterations of the 
beam level optimization during cycle 
No. 8. 
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CYCLES 

Fig. 13 Comparison of the cumulative 

constraint of beam 1 as predicted 
by linear extrapolation at the 
system level and calculated in 
full analysis, (a) infeasible 
design start. - (b) feasible design 



Fig. A1 An example of detailed design 

variables y of a stiffened panel. 
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